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Abstract 

An open problem arising in the statistical description of turbulence is related to the theoretical 
prediction based on first principles of the so-called multi-point velocity probability density functions 
(PDFs) characterizing a Navier-Stokes fluid. 

In this paper it will be shown that - based on a suitable axiomatic approach - a solution to 
this problem can actually be achieved based on the so-called inverse kinetic theory (IKT), recently 
developed for incompressible fluids. More precisely, we intend to show, based on the requirement 
that the Boltzmann- Shannon entropy for the s-point velocity PDF (f s ) is independent of the order 
s and is also maximal at all times, that all multi-point PDFs are necessarily factorized in terms 
of the corresponding 1-point velocity PDF (/i). As a consequence the multi-point PDFs usually 
considered for the phenomenological description of turbulence can be theoretically predicted based 
on the knowledge of f\ achieved by means of IKT. 

PACS numbers: 05.20Jj,05.20.Dd,05.70.-a 



I. INTRODUCTION 



In the context of the statistical description of fluids, the problem of the determination of 
multi-point PDFs arises (at least) in two circumstances: 

• the first one occurs in the phenomenological description of turbulence (see for exam- 
ple Monin and Yaglom [1] 1975 and Pope, 2000 2|]). In such a context, in fact, the 
statistical behavior of fluids is often described in terms of statistical frequencies de- 
fined for multi-point velocity spatial increments (however, similar frequencies can be 
established also for other fluid fields, such as vorticity, scalar pressure, temperature, 
etc.). 

• the second one is the so-called Monin-Lundgren hierarchy 0, 4], based on the construc- 
tion of an infinite set of equations for suitable ensemble-averaged multi-point PDFs 
(ML approach). Such a theory should provide, in principle, also a theoretical model 
for the phenomenological description of turbulence and as a consequence be able to 
predict also the precise form of the velocity-difference PDF observed experimentally in 
HIST (homogenous, isotropic and stationary turbulence). The goal the ML approach 
is actually to predict the time evolution of the ensemble average of the 1-point PDF, to 
be defined in terms of a suitable (and yet to be defined) ensemble-averaging operator. 

Several open issues are related to the ML approach. These concern, in particular, the 
search of possible exact particular solutions of the ML hierarchy represented by a finite set 
of multi-point PDFs. It is well known that the construction of " closure conditions" of this 
type for the ML hierarchy (closure problem) remains one of the major unsolved theoretical 
problems in fluid dynamics. In practice, however, the program of constructing (exact) 
theories of this type or (in some sense) approximate, and holding for arbitrary fluid fields, is 
still open due to the difficulty of preserving the full consistency with the fluid equations. In 
fact, it is well known that many of the customary statistical models adopted in turbulence 
theory - which are based on closure conditions of various type - typically reproduce at most 
only in some approximate (i.e., asymptotic) sense the fluid equations. 

This leaves fundamentally unsolved the problem of the construction of a consistent the- 
oretical model for the multi-point PDFs arising in the phenomenological description of tur- 
bulence. 
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The goal of this paper is to prove (see THM.l below in Section 2) that under suitable 
assumptions all multi-point velocity PDFs characterizing a turbulent NS fluid are factorizable 
in terms of the corresponding 1 -point velocity PDF. 

As a result (see Sec. 3) the treatment of multi-point PDFs can be reached in the con- 
text of IKT (inverse kinetic theory |s 9 ] ) based on the 1-point velocity statistics. It follows 
the fundamental consequence that the multi-point PDFs usually considered for the phe- 
nomenological description of turbulence can actually be theoretically predicted in this way! 
In particular, in the case of local Gaussian 1-point PDF 10] this permits to achieve ex- 
plicit analytic representations of the multipoint velocity PDFs usually considered in the 
phenomenological description of turbulence. 



II. MULTI-POINT STATISTICAL MODEL 



The description of fluids, and more generally of continua, is based on the introduction of 
a suitable set of fluid fields {Z} = {Z iy i — 1, k} satisfying a closed set of PDEs denoted as 
fluid equations. In the case of a fluid obeying of the incompressible Navier-Stokes equations 
(INSE, NS fluid), they are {Z} = {p ,Y,p 1 , St} ■ In particular, here p (the mass density) 
and St (the thermodynamic entropy) are both assumed constant in Q x J, where the latter 
requirement implies that for isentropic flows the equation dST{t)/dt = must hold identi- 
cally for all t G I . In addition V and p\ denote respectively the fluid velocity and the kinetic 
pressure; in particular, pi is defined as the strictly positive function 

Pl {r,t)=p{r,t)+ Po {t) + <f>{r,t), (1) 

where p(r, t),p (t) and 0(r, t) represent respectively the fluid pressure, the (strictly-positive) 
pseudo-pressure and the (possible) potential associated to the conservative volume force 
density acting on the fluid [see the Appendix, Eq.( pH]) ]. 

The statistical description usually adopted for turbulent flows consists, instead, in the 
introduction of appropriate axiomatic approaches denoted statistical models, i.e., sets {/, T} 
formed by a suitable probability density function (PDF) and a phase-space T (subset of 
R n ) on which / is defined. By definition, a statistical model {/, T} realizes a statistical 
description of the fluid if it is possible to define a mapping 

{f,T}^{Z}, (2) 
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which allows the representation in terms of / either: 

A) of the complete set or more generally only B) of a subset of the fluid fields {Z} = 
{Zi,i = l,n} which define the fluid state. 

In particular, the fluid fields Zi(r,t) G {Z} are assumed as functionals of / represented 
by suitable "velocity" moments (of /). In both cases their construction involves, be- 
sides the specification of the phase space (T) and the probability density function (PDF) 
/, the identification of the functional class to which / must belong, denoted as {/} . Sta- 
tistical approaches fulfilling either property A or B will be denoted respectively complete 
and incomplete statistical models. For definiteness in the remainder we shall consider only 
complete statistical models. 

'A priori' the PDF / to be used in a statistical model of this type may be identified with 
an Appoint PDF of the form 

f N (x,t) = /aKxi,...^,*), (3) 
and required to satisfy the normalization condition: 

/ II d\ 3 f N ( X ,t) = l, (4) 

i.e., to be a velocity probability density (in the velocity space U N ); moreover, N > 1 and 
for alH = 1, N, Xj = (r^, Vj) , and Vj denote A" position and velocity vectors, respectively 
belonging to the configuration space of the fluid Q and a suitable velocity space U to be 
identified with M 3 . In particular, consistent with the physical requirement of a NS fluid [i.e. 
the existence of a strong solution of INSE in the set fix/], the following assumptions are 
introduced for f = f N : 

• Axiom #1 (symmetry condition): /jv(xi, x^, t) is symmetric w.r. to arbitrary 
permutation (x 1; ...,Xjy) , i.e., satisfying the invariance condition 

/iv((xi, Xjy) , t) = M(xi, xjv)' , t); (5) 

• Axiom #2 (reduced s-body PDFs): /at(xi, Xjy, t) defines for all s = 1,JV — 1 the 
reduced s-body PDFs 

/ s (xi,...,x s ,t) = — J d 3 r s+1 J d 3 v s+1 / a+ i(xi,...,x s+1 ,i), (6) 

n u 
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where = / <i 3 ri is assumed finite and > 0. Hence each f s satisfies, thanks to @, 
u 

the normalization 

J Hd 3 ^ f s (x 1 ,...,x s ,t) = l; (7) 

Axiom #3 {fluid moments): Jn( x i, •••> x at, determines uniquely the local fluid fields. 
Thus, introducing suitable weight functions G^T}., Vfc,t), for all k = l,n the local fluid 
fields Zi(r k ,t) to be identified with V,pi [both evaluated at the local position and 
time t belonging to Q x I] are taken of the form: 



u (nw-i In?*- 1 11 d 3 r h j uN f] d 3 VjGi(r k , v k ,t)f N (x, i) — 

h=l,N;h^k j=i,N (8) 

= t?Vj/7i(r fc ,v fc ,t)/i(r fc ,v fc ,t) = Zi(r k ,t). 

As suggested by classical statistical mechanics (CSM) [ill ll^j]. Gj(r, v,t) are identified 
respectively with 

G f (r,v,t) = v,p n 2 /3 (9) 
[with u = v — V(r,t) the relative velocity] for V(r,i) and pi(r,i); 

Axiom #4 (entropy moments): ,/jv(t) = /tv( x 15 •••> x at, determines uniquely the 
global fluid field Sr(t). Again based on CSM, the thermodynamic entropy Sr(i) can 
be identified with the Boltzmann-Shannon (BS) statistical entropy. For this reason, 
consistent with Ref. [zj we require that for all t G I: 

Sr(f) = S(/i(f)), (10) 

where /i(t) = /i(xj,t) and /^(x^t) is defined in terms of /jv(£) by means of Eq.(jSJ). 
Furthermore we impose also that for arbitrary iV G Ni and t £ I 

K 2 N S(f N ) = S(h) (11) 

(entropy constraint). Here, denoting by T N the product phase-space T N = [] Ti, 

i=l,N 

with Fx = Q X U, the BS entropy for the iV-point PDF f N is defined as 

S(f N ) = - f dxf N \nf N , (12) 

where dx = FJ dr k dv k and are suitable constants independent /jy to be 

fc=l,iV 

determined; 



• Axiom ^5 (entropic principle): for all iV G Ni, /at(xi, ...,x N ,t) satisfies the principle 
of entropy maximization requiring 

5S(h) = (13) 

(PEM variational principle [l3]]). The variational principle (TIB"]) is imposed either 
solely subject to Axiom #5a (local entropic principle) at some initial time t = t Q or 
to Axiom #5b (global entropic principle) for all t £ I. 

Let us analyze the physical interpretation of the previous assumptions. 

First we notice that #1,#2,#3 and #4 follow from the requirement that the state of the 
fluid is solely prescribed by the set of local and global fluid fields {Z}. In particular the 
locality of the fluid fields, together with the assumption that they are defined everywhere 
in Q, implies manifestly the symmetry requirement (j5]) (see #1). In fact, the positions 
r 1; r N can be manifestly interchanged arbitrarily among them (and similarly the velocity 
vectors v 1; Vjy) without affecting the determination of the fluid fields. This justifies the 
definitions given above in terms of the 1-point PDF both for the local and global fluid fields 
[see Eqs.flED, (DED and TO- 

In a similar way, since by assumption the fluid fields cannot depend on the level adopted 
for the statistical description of the fluid, all moments which define the fluid field must be 
independent of the choice of the iV-point PDF. Indeed, for arbitrary N £ Nj., it must be 
possible to represent the thermodynamic entropy in terms of Boltzmann- Shannon entropy 
associated to the iV-point PDF f^, as well to f\. This implies, that besides the position 
(jHJ) invoked for the local fluid fields also the additional constraint ffTTl) must be placed on 
all the BS entropies associated to multi-point PDFs. This constraint manifestly should hold 
identically (for all t £ J). 

Finally, the hypothesis that the Boltzmann-Shannon entropy is maximal (see #5) implies 
the validity of the entropic principle (fTBl) . We stress that, in principle, PEM can be assumed 
to hold either at the initial time t or, more generally, for arbitrary t £ /. The second 
requirement is consistent with the assumption of isentropic flow. In fact, the positions (flQj) 
and fflTl) imply that also the BS entropy must be constant (i.e., independent of time). Hence, 
the requirement that it is maximal at some initial time t may not be at variance with the 
requirement placed by the global entropic principle #5b. 
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Basic issues are related to the, possibly non-unique, determination of the appropriate 
statistical model {/, T}. These concern in particular: 

1. (PROBLEM j^l) the search of the (possible) minimum level (N) of the statistical 
description to be adopted for {/, T} ; 

2. (PROBLEM #2) the determination of the time-evolution of the multi-point PDFs fx] 

3. (PROBLEM #3) the determination of the initial and boundary conditions for f^. 

Regarding the first problem the following remarkable result holds: 

THM.l - Factorization theorem for f N . 

Let us impose Axioms #l-#4 with #5b. Then it follows necessarily that: 

1) the variational constraint 

5{K 2 N S(f N )-S(f 1 )} = (14) 

must hold for all t £ I; 

2) for all N G Ni, the iV-point PDF f N (xi, ...,x N ,t) is of the form: 

fa(x 1 ,...,x N ,t)= Yl /i( x i>*)> ( 15 ) 

i=l,N 

with /i(xi,£) denoting the corresponding 1-point PDF defined by Eq.§B\). Hence, it follows 
also that for all s — 1, N — 1; 

f s (x 1 ,...,x s ,t) = J d 3 v s+1 / s+1 (x!,...,x s+1 ,t). (16) 

u 

3) the constant in Eg. /[Jl\) reads 

K 2 N = Nfi(nf-\ (17) 

PROOF First we notice that the entropy constraint (1111) together the global entropic prin- 
ciple #5b [i.e., the requirement that Eq. ffTB"]) holds for all t G I] imply that, for all iV 
and for all t G /, also the variational constraint (Tl4"|) must be fulfilled. To prove that the 
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factorization property of the iV-point PDF must hold for all t E I, let us consider for illus- 
tration (and without loss of generality) the case N = 2. Denoting / 2 (xi, x 2 , i) = / 2 (1, 2) and 
/i(xi,£) = /i(l), Eq. fH]) delivers for arbitrary variations Sfi(3): 

J dx5A(3) {/ 2 (l,2)ln/ 2 (l,2) -A(l)/i(2) [In h (1) + In A (2)]} = 0. (18) 
r 3 

This implies necessarily that the factorization condition / 2 (1,2) = / 1 (1)/ 1 (2) must hold 
identically in V 2 x /. The proof can easily be extended to arbitrary N > 2, yielding Eq-f TTS]) . 
In turn, thanks to Eq. f!15j) . equations f )16p and (\17\i immediately follow, respectively from 
Eqs.flHD and (JIT). Q.E.D. 

We remark that in principle THM.l can be generalized by requiring that PEM holds only 
at the initial time t Q G / (Axiom #5a). Nevertheless, in this case the constraint ( TTT|) only 
warrants that the factorization condition ffTB 7 ]) holds at the initial time t Q , unless the form 
of the statistical (Liouville) equations holding for the s-point velocity PDFs is explicitly 
prescribed as done in Ref. |9(. 

Invoking, however, the validity of Axiom #5b and consequently of THM.l, the statistical 
model {/, T} can be identified with the IKT statistical model for the 1-point PDF 



III. IKT FOR MULTI-POINT PDFS 



The construction of multi-point PDFs is a problem of "practical" interest in experi- 
mental/numerical research in fluid dynamics, usually adopted for the statistical analysis of 
turbulent fluids. In fact, they can be experimentally measured in terms of velocity differences 
between different fluid elements. 

Let us assume, for defmiteness, that /i(xj,t) is the 1— point PDF which is particular 
solution of the Liouville equation [or inverse kinetic equation (IKE)] provided by IKT js|. 
Then, denoting f\(i) = A(xj,t) (for i = l,s) the same PDF evaluated at the states Xj = 
(rj, Vj) (for i — 1, s), the s— point PDF is the probability density 

/ S (1,2,.. S )= J] (19) 

1=1, s 

defined in the product phase-space V s = ]~[ T, The statistical equation advancing in time 

i=l,s 

f s follows trivially from the Liouville equation for the 1-point PDF see [5]). In fact, denoting 
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for i = 1, s by F(i) = F(xj, fi) the 1-point mean-field force per unit mass acting on the i-th 
particle (with state Xj) [defined in Refs. Q and Q] and introducing the s— point Liouville 
operator 



L s (l,..,s 



8=1, S 



<9r, : <9r, ; 



(20) 



it follows that / s (l,2, ..s) satisfies identically the s— point Liouville equation 



L a (l,..,s)/,(l,2,..s) = 0. 



(21) 



A. Explicit evaluation of 2-point velocity PDFs 

In terms of the 2-point PDF, /2(1,2), a number of reduced probability densities can 
be defined in suitable subspaces of T 2 . To introduce them explicitly let us first introduce 
the transformation to the center of mass coordinates of the two point-particles with states 
(T t , V;) (for % = 1,2) 

{ ri , Vl ,r 2 ,v 2 }^{r,R,v,V} (22) 

[here r = , R = ri ^ r2 ; furthermore, v, V can be identified with v = vi — V2 and 
V = v x + v 2 ]. Then, these are respectively: 

1) the local (in configuration space) velocity-difference 2-point PDF #2(1*1, r 2 , v, t) defined 
in the phase-space Q 2 x U and obtained integrating the 2-point velocity PDF w.r. to the 
mean velocity V 

fla(r lj ri l v l t) = / E ,d 8 V/ 2 (l l 2) = 
= / d 3 VA( ri , v + V,t))A(r 2 , V - v,t); 

2) the velocity-difference 2-point PDF / 2 (r,v,t) defined in Ti = Q X U and obtained 
integrating also on the center-of-mass position vector R. Thus denoting by 

1 r j3 , 



<-) a , n = ^ I (24) 



n 



the configuration-space average operator acting on the center of mass coordinates R, 
there it follows 

/ 2 (r, v,t) = (g 2 (r + R, R - r, v, t)) R Q . (25) 



9 



In particular, in the case of a Gaussian PDF 9j, Eq. (l23|) delivers again a Gaussian-type 
PDF 



2 (ri,r 2 ,v,t) 



1 



n 3 / 2 v'f h 



exp 



V(l)-V(2) 



(26) 



■tii 



where V(z) = V(rj,t), v. 



th,p\ 



J th,p 



(rj,t) and v\ h denotes 



<p(l) + < P ( 2 ) 
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(27) 



In a similar way it is possible to obtain explicit representations for the following additional 
2-point PDFs: 

1. the velocity- difference 2-point PDF for parallel velocity increments . Introducing the 
representations v = nu and r = nr, n denoting a unit vector, f 2 u(r,v,t) can be simply 
defined as the solid-angle average 



f 2 \\{r,v,t) = / dCl(n)f 2 (r = nr, v = nv,t); 



(28) 



2. the velocity- difference 2-point PDF for perpendicular velocity increments . Introducing, 
instead, the representations v = nv and r = n x br, n and b denoting two independent 
unit vectors, f2±(r,v,t) can be defined as the double-solid-angle average 

(29) 



f2±(r,v,t) = J dfi(n) J dfi(b) 
f 2 (r = n x br, v = nv,t). 

An interesting property which emerges from these results is that in all cases indi- 
cated above [i.e., Eqs.( l25]) .( l28]) and (|29|) ] the definition of g 2 given above [Eq. ( 123]) ] 
implies that non-Gaussian features, respectively in f 2 ,f 2 \\ and f 2 ±, may arise even if 
the 1— point PDF is Gaussian. This occurs due to velocity and pressure fluctuations 
occurring between different spatial positions ri and r 2 . More generally, however, we 
can infer that, due to the constraint here imposed on the 1-point PDF 



(h{t)) r n = W req \t) 



(30) 



[where (-) rn it the averaging operator (•) rQ = -^y / d 3 r acting on of a function 

n 

F(x, t)}, it is obvious that, if the fluid velocity V(r, t) is bounded in the domain Q, the 
same 1-point PDF, and hence the 2-point PDFs, cannot be Gaussian distributions. 
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B. Statistical evolution equation for the velocity-difference 2-point PDF 



From the 2— point IKE (12 ip (obtained in the case s = 2) it is immediate to obtain the 

corresponding evolution equation for the reduced PDFs indicated above. For example, the 
velocity-difference 2-point PDF f 2 satisfies the equation 

d h , 9 ? 9 r> 

^ + V ^ /2 = -^' D (31) 

where D is the diffusion vector 

D= f <PV ( Fl{1) ~ F2{2) f 2 (l,2)) . (32) 



It follows, in particular, that in the case of a Gaussian 1-point PDF this equation reduces 
to the Fokker-Planck equation 

9 fo | v 9 ? = 9 p (33) 
dt dr dv 

where the Fokker-Planck diffusion vector D reads 

D = (F^ 2 (r + R, R - r, v, *)) Rn (34) 

and the vector field = F^(ri, r 2 , V, I-Jm) is reported in Ref. 0. It follows that both 
equations are manifestly non-Markovian as a consequence of the non-local dependencies 
arising (in both cases) in the Fokker-Planck coefficients D and D. 



An interesting issue is here provided 
veloped by Peinke and coworkers 
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)y the comparison with the statistical formulation de- 
18| . Their approach, based on the statistical analysis of 
experimental observations, indicates that in case of stationary and homogeneous turbulence 
both the 2-point PDFs for parallel and velocity increments obe y st ationary Fokker-Planck 
equations. In particular, according to experimental evidence 13, |l8| a reasonable agreement 
with a Markovian approximation for Eq. (I33j) - at least in some limited subset of parameter 
space- is suggested. Our theory implies, however, that a breakdown of the Markovian prop- 
erty should be expected due to non-local contributions appearing in the previous statistical 
equations (|3Tjl and (|33]). 

IV. CONCLUSIONS 

In this paper we have shown that the multi-point PDFs used in customary phenomeno- 
logical approaches to turbulence can be explicitly evaluated in terms of the 1-point velocity 
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PDF (fx) determined in the framework on the IKT-statistical model 0-9]. 

The starting point is provided by THM.l, which shows that under suitable hypotheses 
the multi-point PDF is necessarily factorized in terms of the 1-point PDF fx. The re- 
quirements here imposed include, in particular, the assumption that {/at, T n } is a complete 
statistical model, i.e., that in terms of the multi-point PDF the complete set of fluid fields 
(defining the fluid state) can be represented by means of suitable velocity and phase-space 
moments [see Axioms $4-^6]. Then, provided: 

A) the entropy constraint (ITT|) is invoked ((Axiom #4); 

B) the validity of PEM is imposed at all times t 6 I (Axiom #5b); 

the factorization condition (II 5p for in terms of the 1-point PDF fx necessarily follows. 

As a result, in validity of the previous requirements, the statistical model for NS fluid 
can be identified with the IKT-statistical model {/i,ri} earlier developed [5 - 8 ] and based 
on the 1-point PDF fx- The theory has important consequences: 

1. arbitrary multi-point PDFs can be uniquely represented in terms of the 1-point PDF 
characterizing the IKT-statistical model {/i,ri}; 

2. the time evolution of the multi-point PDFs is uniquely determined by {fx, Tx} ; 

3. the theoretical prediction of multipoint PDFs is actually possible. 

4. qualitative properties of the multi-point PDFs can be investigated. As a particular 
case, the example of a Gaussian 1-point PDF has been pointed out. 

In the IKT-statistical model the statistical equation advancing in time the 1— point PDF 
fx coincides with the Liouville equation. As a consequence, its explicit evaluation is actually 



made possible js]. In particular, as shown in Ref. [lOj, in the presence of HIST the 1-point 
PDF necessarily coincides with a Gaussian distribution. Thanks to the factorization theorem 
(THM.l) this implies that also the multi-point velocity PDFs are uniquely determined. 
As result, as indicated in Section 3 (see subsection 3.1), two-point PDFs relevant for the 
phenomenological description of hydrodynamic turbulence can be explicitly determined. 
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V. APPENDIX: INSE PROBLEM 

The fluid equations for a NS fluid are the so-called incompressible Navier-Stokes equations 
(INSE) for the fluid fields {Z} = {p , V(r,t),p x (r,t), S T (t)} : 



P = Po, (35) 

V- V = 0, (36) 

NV = 0, (37) 

= 0, (38) 

Z{r,t ) = Z {r), (39) 

Z(r,t)\ m = Z w (r,t)\ m , (40) 



Eqs. ( 1351) - f j40l) denote respectively the incompressibility, isochoricity, Navier-Stokes and 
constant thermodynamic entropy equations and the initial and Dirichlet boundary conditions 
for {Z} , with {Z a (r)} and {Z w (r,t)\ gn } suitably prescribed initial and boundary-value fluid 
fields, defined respectively at the initial time t = t Q and on the boundary dfl, In particular, 
this means that they are are required to be at least continuous in all points of the closed set 
fl x I, with fl = fl U dQ closure of Q. In the remainder we shall require that: 

1. Q (configuration domain) is a bounded subset of the Euclidean space E 3 on M 3 ; 

2. I (time axis) is identified, when appropriate, either with a bounded interval, i.e., 
I=]to, tx[ C R, or with the real axis R; 

3. in the open set Qxl the functions {Z} , are assumed to be solutions of Eqs. fl36|) - fl38|) 
subject, while in fix I they satisfy the whole set of Eqs. (!35|) - P0|) . In particular: Eqs. 
( |35l) - (HOJ) define the initial-boundary value INSE problem, 

4. by assumption, the fluid fields are strong solutions of the fluid equations. Hence Eqs. 
fl35|) - f|40|) are required to define a well-posed problem with unique strong solution 
defined everywhere in fix I. 

Here the notation as follows. iV is the NS nonlinear operator 

NV=^-V-F H , (41) 
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with jy: V and F# denoting respectively the Lagrangian fluid acceleration and the total force 
per unit mass 

= ^V(r,t) + V(r,t) ■ VV(r, t), (42) 
F H = --Vp(r,t) + -f (r,t) + t;V 2 V(r,t), (43) 

while p > and v > are the constant mass density and the constant kinematic viscosity. 
In particular, f is the volume force density acting on the fluid, namely which is assumed of 
the form 

f = -V0((r,t) + f fl (r,t), (44) 

<f)((r,t) being a suitable scalar potential, so that the first two force terms [in Eq. (]43l ] can be 
represented as 

- Vp(r,t) + f (r,t) = -V Pl (r,t) + f R (r,t), (45) 

with pi(r,t) defined by Eq.ffl]) denoting the kinetic pressure. As a consequence the fluid 
pressure necessarily satisfies the Poisson equation 

V 2 p(r,t) = S(r,t), (46) 

where the source term S reads 

S(r,t) = -p V • (V • VV) + V • f . (47) 
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